Black-hole quasinormal resonances: Wave analysis versus a geometric-optics 

approximation 



Shahar Hod 

The Ruppin Academic Center, Emeq Refer 40250, Israel 

and 

The Hadassah Institute, Jerusalem 91010, Israel 
(Dated: September 2, 2009) 

It has long been known that null unstable geodesies are related to the characteristic modes of 
black holes- the so called quasinormal resonances. The basic idea is to interpret the free oscillations 
of a black hole in the eikonal limit in terms of null particles trapped at the unstable circular orbit 
and slowly leaking out. The real part of the complex quasinormal resonances is related to the 
angular velocity at the unstable null geodesic. The imaginary part of the resonances is related to 
the instability timescale (or the inverse Lyapunov exponent) of the orbit. While this geometric- 
optics description of the black-hole quasinormal resonances in terms of perturbed null rays is very 
appealing and intuitive, it is still highly important to verify the validity of this approach by directly 
, analyzing the Teukolsky wave equation which governs the dynamics of perturbation waves in the 

black-hole spacetime. This is the main goal of the present paper. We first use the geometric- 
optics technique of perturbing a bundle of unstable null rays to calculate the resonances of near- 
extremal Kerr black holes in the eikonal approximation. We then directly solve the Teukolsky 
wave equation (supplemented by the appropriate physical boundary conditions) and show that the 
■ resultant quasinormal spectrum obtained directly from the wave analysis is in accord with the 

spectrum obtained from the geometric-optics approximation of perturbed null rays. 
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The no-hair conjecture [l| asserts that the external held of a perturbed black hole relaxes to a Kerr-Newman 
$H , spacetime, characterized solely by three parameters: the black-hole mass, charge, and angular momentum. This 
^Q 1 implies that perturbation helds left outside the black hole would either be radiated away to infinity, or be swallowed 
by the black hole. 

The relaxation phase in the dynamics of perturbed black holes is characterized by 'quasinormal ringing', damped 
oscillations with a discrete spectrum (see e.g. 0] for reviews). At late times, all perturbations are radiated away in 
a manner reminiscent of the last pure dying tones of a ringing bell 0, 0i S 0| ■ Quasinormal resonances are expected 
\ to play a prominent role in gravitational radiation emitted by a variety of astrophysical scenarios involving black 
holes. Being the characteristic 'sound' of the black hole itself, these free oscillations are of great importance from the 
\ astrophysical point of view. They allow a direct way of identifying the spacetime parameters, especially the mass and 
■ angular momentum of the black hole. 

The dynamics of black-hole perturbations is governed by the Regge- Wheeler equation Q in the case of a spherically 
symmetric Schwarzschild black hole, and by the Teukolsky equation Q for rotating Kerr-Newman spacetimes. The 
\ black hole quasinormal modes (QNMs) correspond to solutions of the wave equations with the physical boundary 
conditions of purely outgoing waves at spatial infinity and purely ingoing waves crossing the event horizon |9| . Such 
boundary conditions single out a discrete set of black- hole resonances {to n } (assuming a time dependence of the form 
e~ tuJt ). In analogy with standard scattering theory, the QNMs can be regarded as the scattering resonances of the 
black-hole spacetime. They thus correspond to poles of the transmission and reflection amplitudes of a standard 
scattering problem in a black-hole spacetime. 

In accord with the spirit of the no-hair conjecture the external perturbation fields would either fall into the 
black hole or radiate to infinity. This implies that the perturbation decays with time and the corresponding QNM 
frequencies are therefore complex. It turns out that there exist an infinite number of quasinormal modes, characterizing 
oscillations with decreasing relaxation times (increasing imaginary part), see 13, he mm 

and references therein. 
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The mode with the smallest imaginary part (known as the fundamental mode) determines the characteristic dynamical 
timescale for generic perturbations to decay [3 El, HI, 17, H, ljj ]. 



In most cases of physical interest [20( , the black-hole QNMs must be computed numerically by solving the Teukolsky 
equation supplemented by the appropriate physical boundary conditions. However, Mashhoon [2l| has suggested an 
analytical technique of calculating the QNMs in the geometric-optics (eikonal) limit. The basic idea is to interpret 
the black-hole free oscillations in terms of null particles trapped at the unstable circular orbit and slowly leaking out 
[2H I22I [23j . The real part of the complex quasinormal resonances is related to the angular velocity at the unstable 
null geodesic while the imaginary part of the resonances is related to the instability timescale of the orbit (or the 
inverse Lyapunov exponent of the geodesic [23|). 

It should be emphasized that, Mashhoon's approach of computing black- hole quasinormal frequencies in the eikonal 



2 



limit is somewhat indirect: instead of explicitly solving the Teukolsky wave equation which describes the dynamics 
of perturbation waves in the black-hole spacetime, Mashhoon analyzed perturbations of test null rays in the unstable 
circular orbit of the black hole. While this geometric-optics description of black-hole quasinormal resonances in terms 
of perturbed null rays is very appealing and intuitive, it is still important to verify the validity of this approach 
by directly solving the Teukolsky master equation which governs the dynamics of perturbation waves in the Kerr 
black-hole spacetime. Filling this gap is the main goal of the present paper. We shall first use the geometric-optics 
technique of perturbing a bundle of unstable null geodesies to calculate the resonances of near-extremal Kerr black 
holes in the eikonal limit. We shall then directly solve the Teukolsky wave equation (supplemented by the appropriate 
physical boundary conditions) and show that the resultant resonances of the wave equation agree with those obtained 
from Mashhoon's indirect approach of perturbing null rays in the appropriate eikonal approximation. 

We start with Mashhoon's approach of calculating the black- hole QNM resonances in the eikonal limit I = m ^> 1, 
where I is the angular-momentum of the perturbed null rays, and m is the azimuthal harmonic index of the rays. 
Accordin g to Mashhoon's analysis of perturbed null geodesies, the Kerr QNM frequencies in the I — m ^> 1 limit are 
given by [2l| (we use natural units in which G = c = % = 1) 

u n = mu)+ - i(n+^)0u) + ; n = 0,1,2,..., (1) 

where 

(2) 



r%* + aMl/2 



is the Kepler frequency for null rays in the unstable equatorial circular orbit of the black hole, and 

r ph = 2M{1 + cos[^ cob-^-o/M)]} , (3) 



is the limiting circular photon orbit. The function (3 is given by 



(12M) 1 / 2 (r p /i — r + )(r ph — r_) 



p= K - — ' lii: + \\r . ( 4 ) 

r ph 



rT(r ph -M) 



see Ref. |2l| for details. Here M and a are the mass and angular momentum per unit mass of the black hole, 
respectively. 

We now focus on the near extremal limit Tbh ^ 0, where 

(M 2 - a 2 ) 1 / 2 f . 

BH AttM[M + (M 2 — a 2 ) 1 / 2 ] ' 1 j 

is the Bekenstein-Hawking temperature of the black hole. Let r± — M±e with e/M <C 1, where r± = M +(M 2 — a 2 ) 1 ' 2 
are the black- hole event and inner horizons. The Bekenstein-Hawking temperature of the black hole now reads 
Tbh = e/4?rAf 2 + 0(e 2 /M 3 ). After some algebra we find from Eq. © 



This, in turn, implies 



and 



r ph = M+^= + 0{e 2 /M) . (6) 



u + = — -^ + 0(e 2 /M 3 ) , (7) 
+ 2M AM 2 y ' ' ' V ; 



P=- + 0(e 2 /M 2 ). (8) 

Substituting Eqs. (JT]) and © into Eq. (fT]), one finds that the black- hole quasinormal resonances in the eikonal 
approximation I = m >> 1 are given by 

fz 1 

uj n = mn + 2ttT B h(1 - \)m - i2^T BH {n + -) + 0(MT 2 H ) , (9) 



3 



where 

n = t\t^ = wr 2nTBH + °( MT ^ ' ( 10 ) 

is the angular velocity of the black-hole event horizon. 

We shall now study directly the Teukolsky wave equation in order to determine the fundamental (least-damped) 
resonant frequencies of the Kerr black hole. The Teukolsky equation is amenable to an analytical treatment in the 
near-extremal limit (M 2 - a 2 ) 1 /' 2 <C a < M. 

In order to determine the black-hole resonances we shall analyze the scattering of massless waves in the Kerr 
spacetime. The dynamics of a perturbation field \f in the rotating Kerr spacetime is governed by the Teukolsky 
equation [8J. One may decompose the field as 

¥ rim (t, r, 9, 4,) = e im<t, S s i m (6; auj)^ slm (r)e- wt , (11) 

where (t,r,9,<f>) are the Boyer-Lindquist coordinates, u> is the (conserved) frequency of the mode, / is the spheroidal 
harmonic index, and m is the azimuthal harmonic index with —l<m<l. The parameter s is called the spin weight 
of the field, and is given by s = ±2 for gravitational perturbations, s = ±1 for electromagnetic perturbations, s = ±A 
for massless neutrino perturbations, and s = for scalar perturbations. (We shall henceforth omit the indices s,l,m 
for brevity.) With the decomposition (jll[) . ip and S obey radial and angular equations, both of confluent Heun type 
[U^, coupled by a separation constant A(au). 

The angular functions S(9; auj) are the spin-weighted spheroidal harmonics which are solutions of the angular 
equation pi EHl] 



1 dr. BS^ 



(m + s cos ( 
sin 2 



( sin 9^-) + \a 2 u 2 cos 2 9 - 2au JS cos 9 - l m + 8C0Sf/ J + s + A ] g = . (12 ) 
>d9\ 89 J L sin 2 9 1 



The angular functions are required to be regular at the poles 9 — and 9 = ir. These boundary conditions pick out a 
discrete set of eigenvalues Ai labeled by an integer I. [In the au> 1 limit these angular functions become the familiar 
spin-weighted spherical harmonics with the corresponding angular eigenvalues A — 1(1 + 1) — s(s + 1) + O(atu).] The 
angular equation (IT2|) can be solved analytically in the I = m >> 1 limit to yield 



A = m 2 +0(m) , (13) 



seeRef. (26j. 

The radial Teukolsky equation is given by 



. , d/ A ,4.i#\ \K 2 — 2is(r — M)K 2 2 „ A 

A- s — (A s+1 -^) + \ a 2 uj 2 + 2mau; - A + Aiswr 

dr\dr/lA 



rl> = 0, (14) 



- 2Mr + a 2 and K = (r 2 +a 2 )uu — am. For the scattering problem one should impose physical boundary 
conditions of purely ingoing waves at the black-hole horizon and a mixture of both ingoing and outgoing waves at 
infinity (these correspond to incident and scattered waves, respectively). That is, 

( e -^v + K{uj)e iuJV as r -> 00 (y -> 00) ; 
~ \T(oj)e- l ^~ m ^ y as r r+ (y -> -00) , ^ ' 

where the "tortoise" radial coordinate y is defined by dy = \(r 2 + a 2 )/A]dr. The coefficients T{uS) and lZ(u>) are 
the transmission and reflection amplitudes for a wave incident from infinity. They satisfy the usual probability 
conservation equation |T(w)| 2 + |7?.(w)| 2 = 1. 

The discrete quasinormal frequencies are the scattering resonances of the black-hole spacetime. They thus corre- 
spond to poles of the transmission and reflection amplitudes. (The pole structure reflects the fact that the QNMs 
correspond to purely outgoing waves at spatial infinity.) These resonances determine the ringdown response of a black 
hole to external perturbations. Teukolsky and Press [13] and also Starobinsky and Churilov [111 have analyzed the 
black-hole scattering problem in the double limit a M and u> — * mfL Detweiler [2j| then used that solution to 
obtain a resonance condition for near-extremal Kerr black holes. Define 

T-t- — T — 

a = ; t = Mho — mVl) ; <1> = ujr + . (16) 
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Then the resonance condition obtained in [29[ for a << 1 and r << 1 is: 

T(2i6)T(l + 2i6)T(l/2 + s-2iu-id)T(l/2-s-2iu-i6) _ ^ 2i ~.^2i5 r (V 2 + 2 ^ + iS ~ ^r/a) ^ 



T(-2iS)T(l-2iS)T(l/2 + s-2ii2} + iS)T(l/2- s-2iuj + i5) K ' Y{l/2 + 2iuo - iS ~ Air/a) ' 
where 8 2 = Alj 2 — 1/4 — A — a 2 u> 2 + 2maui. Taking cognizance of Eq. (fT3")) , one finds 

S =^ m + 0(1), (18) 

for near-extremal Kerr black holes in the I = to >> 1 limit. 

The left-hand-side of Eq. (TTT)) has a well defined limit as a — > M and w — ► mQ. We denote that limit by C. Since 
S ~ >> 1, one has (-i)- 2l<5 = e-iVSmM-i) = e -iV3mine-™/ 2 = e -iVSm(-^/2) = g-VSwn/s < 1; whicn 

implies A = (— 2ia)er)~ 21 ' 5 -C 1. Thus, a consistent solution of the resonance condition, Eq. (TT7|) . may be obtained if 
l/T(l/2 + 2iu -iS- Air/a) = 0(A). Suppose 

1/2 + 2i(Jj -id- Air/a = -n + r 1 X + 0(A 2 ) , (19) 

where n > is a non-negative integer and r\ is an unknown constant to be determined below. Then one has 

T(l/2 + 2iu -iS- Air/a) ~ r(-n + rjX) ~ (-ri) _1 r(-n + 1 + rye) ~ • • • ~ [(-l) n n!] _1 r(7?A) , (20) 

where we have used the relation T(z + 1) = zT(z) Next, using the series expansion l/r(z) = J2kLi c k zk with 

ci = 1 [see Eq. (6.1.34) of [30j]]. one obtains 

l/r(l/2 + 2iu> -iS- Air/a) = (-l)"n!r;A + 0(A 2 ) . (21) 

Substituting this into Eq. (TTT]) one finds 77 = £/[(-l)"n!r(-n + 2iS)]. 

Finally, substituting At /a = (u — mf2) /2t:Tbh , 2id) = im + O^tuMTbh) f° r w — + 0(77iTbh), and (5 = 
V3m/2 + O(l) into Eq. ((Til), one obtains the resonance condition 

/3 

(lj -mn)/2'!tT B H = i[-n + r}\- l/2 + «m(-y - 1)] . (22) 
The black- hole quasinormal resonances in the I = to >> 1 limit are therefore given by the formula 

w„ = toO + 27rT BH (l - ^_)m - i2^T BH {n + -) + 0{MT 2 BH ) , (23) 

where n — 0,1,2,... . We emphasize that this result, obtained from the direct wave analysis, coincides with the 
previously derived spectrum Eq. ^ of the ray analysis (the geometric-optics approximation). 

In summary, we have studied analytically the quasinormal spectrum of rapidly-rotating Kerr black holes. We first 
used the technique of perturbing a bundle of unstable equatorial null geodesies to calculate the quasinormal resonances 
of Kerr black holes in the cikonal limit I = m » 1. We then used the alternative (and the more direct) approach 
of solving the Teukolsky wave equation which governs the dynamics of perturbation waves in the Kerr black-hole 
spacetime. We have shown that the resonance spectrum (I23|) obtained directly from the wave analysis is in accord 
with the spectrum ^ which was obtained from the geometric-optics approximation of perturbed null rays. 
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